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Abstract 

We  extend  our  recent  work  where  a  screw  dislocation  in  FCC  Ni  was  found  to 
spontaneously  attain  a  low-energy  partially  cross-slipped  configuration  upon  intersecting 
a  forest  dislocation  [1].  Using  atomistic  (molecular  statics)  simulations  with  embedded 
atom  potentials,  we  evaluate  the  activation  barrier  for  a  dislocation  to  transform  from 
fully  residing  on  the  glide  plane  to  fully  residing  on  the  cross-slip  plane  intersecting  a 
forest  dislocation.  The  activation  energies  were  obtained  by  detennining  equilibrium 
configurations  (energies)  when  variable  pure  tensile  or  compressive  stresses  are  applied 
along  the  [111]  direction  on  the  partially  cross-slipped  state.  We  show  that  the  activation 
energy  is  a  factor  of  3  -  6  lower  than  that  for  cross  slip  in  isolation  via  the  Escaig  process. 
Further,  the  activation  barrier  for  cross-slip  at  these  intersections  is  shown  to  be  linearly 
proportional  to  (d/b)ln(d/b),  as  in  the  Escaig  process,  where  ‘d’  is  the  Shockley  partial 
dislocation  spacing  and  „b’  is  the  Burgers  vector  of  the  screw  dislocation.  These  results 
suggest  that  cross-slip  should  be  preferentially  observed  at  selected  screw  dislocation 
intersections  in  FCC  materials. 


1 


Keywords:  Cross-slip;  Activation  analysis,  Atomistic  simulations;  Nickel;  Escaig 

stresses 

1.0  Introduction 

There  is  an  increasing  recognition  of  the  need  to  incorporate  physics-based  models  of 
deformation  in  design  of  structural  components.  While  models  for  predicting  yield 
strength  and  creep  behavior  are  beginning  to  incorporate  significant  physics,  models  for 
fatigue  and  ultimate  strength  remain  mostly  empirical.  This  is  because  strain-hardening 
and  fatigue  resistance  are  highly  influenced  by  dislocation  micromechanisms,  for 
example  cross-slip,  and  including  physics-based  cross-slip  processes  in  mesoscale 
simulations  (e.g.  dislocation  dynamics)  has  been  difficult.  The  early  work  of  Escaig 
remains  the  most  widely  cited  and  used  model  for  cross-slip  [2-5];  however,  this  model 
poses  several  difficulties  with  respect  to  quantitative  simulations.  The  model  is  highly 
sensitive  to  choice  of  parameters  that  have  thus  far  been  difficult  to  obtain.  For  example, 
the  constriction  width  required  for  cross-slip  is  unknown  but  significantly  influences  the 
enthalpy  of  the  cross-slip  process  [6].  This  difficulty  has  led  to  an  ad-hoc  postulate  that 
obstacles  always  exist  in  materials  which  enable  sufficient  dislocation  core  constriction 
under  stress  thereby  ensuring  cross-slip  [2],  Aside  from  being  a  presumption,  this  forces 
cross-slip  models  (particularly  for  single  crystals)  to  make  arbitrary  assumptions  about 
the  dislocation  obstacle  spacing.  Advances  in  atomistic  simulations  make  it  possible  to 
gain  insights  into  the  cross-slip  process  and  may  serve  to  inform  mesoscale  simulations  to 
accurately  capture  the  atomic-level  physics  of  that  dislocation  process.  Such  an  insight 
was  realized  in  our  recent  work[l],  where  the  possibility  of  cross-slip  at  dislocation- 
dislocation  intersections  was  selectively  examined  using  large-scale  atomistic  simulations 
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that  contain  more  than  one  dislocation.  A  glide  dislocation  that  intersects  a  pair  of  forest 
dislocations  was  modeled  and  several  conditions  were  shown  to  exist  where  cross-slip 
can  nucleate  easily.  This  provides  a  better  physical  basis  for  cross-slip  and  means  to 
incorporate  realistic  cross-slip  evolution  in  larger-scale,  discrete  dislocation-dynamics 
simulations.  Further,  this  may  lead  to  a  realistic  statistical  representation  of  cross-slip 
during  monotonic  or  cyclic  deformation.  However,  the  prior  work  did  not  explore  the 
size  and  nature  of  the  activation  barrier  for  the  dislocation  transformation  at  intersections. 

In  this  work  we  show  that  a  finite  energy  barrier  exists  for  cross  slip  at  an 
intersection,  but  that  it  is  much  smaller  than  that  currently  understood  from  analyses  of 
the  Escaig  process.  Molecular  statics  simulations  were  used  to  determine  the  saddle  point 
configuration  and  activation  energy  by  changing  the  extent  of  partial  dislocation  cross¬ 
slip  using  an  applied  stress.  The  applied  stresses  were  compressive  or  tensile  along  the 
[111]  direction  on  the  partially  cross-slipped  dislocation  structures  obtained  in  prior  work 

[1].  The  screw-character  dislocation  under  study  had  a  lA[  \  1  0]  Burgers  vector  and  the 
intersecting  dislocation  had  the  'A[  101]  Burgers  vector  on  the  (111)  plane  and  a  120° 
line  orientation  of  [0-1-1].  Section  2  describes  the  simulation  technique  and  the 
interatomic  potential  used  for  the  simulations.  Section  3  presents  a  review  of  the  core 
structures  obtained  from  atomistic  simulations  for  this  specific  intersection  [1],  the  effect 
of  applied  stresses  along  the  [111]  direction  on  the  partially  cross-slipped  core  structure, 
the  identification  of  the  saddle  point  configuration  and  a  detennination  of  the  energies  of 
the  activated  state  for  cross-slip.  Finally,  Sections  4  and  5  give  a  discussion  and  summary 
of  the  results  respectively. 
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2.0  Simulation  technique 

The  atomistic  simulations  described  here  employed  the  3 -dimensional  (3d)  parallel 
molecular  dynamics  code,  LAAMPS  [7],  developed  at  Sandia  National  Laboratory.  A 
schematic  of  the  simulation  cell  used  in  the  atomistic  simulations  is  given  in  Figure  1. 
The  simulation  cell  is  a  rectangular  parallelepiped  having  the  x-axis  oriented  along 
[110],  y-axis  along  [112]  and  the  z-axis  along  [111].  The  dimensions  of  the  simulation 
cell  are  62.0  nm  along  the  x-axis  and  31.5  nm  along  both  the  y-  and  z-  axes 
corresponding  to  a  simulation  cell  of  5,405,160  atoms.  The  partially  cross-slipped  core 
structures  obtained  for  the  screw  (l/2[  1  10]  Burgers  vector)  -  120°  (l/2[  101]  Burgers 
vector,  [011]  line  direction)  intersection  in  reference  1  were  subjected  to  unifonn 
compressive  or  tensile  stresses  applied  along  the  [111]  direction.  Under  such  a  state  of 
stress,  there  are  no  resolved  shear  stresses  acting  along  the  Burgers  vector  direction  for 
either  the  screw  or  the  120°  dislocation.  Also,  there  are  no  applied  Escaig  stresses  acting 
on  the  edge  components  of  the  screw-character  dislocation  on  the  (111)  glide  plane. 
However,  there  are  Escaig  stresses  acting  on  the  edge  components  of  the  screw 
dislocation  on  the  cross-slip  (111)  plane  as  well  as  the  forest  dislocation  on  the  (111) 
plane.  Fixed  boundary  conditions  were  applied  along  all  three  directions  and  energy 
minimization  was  perfonned  using  the  conjugate  gradient  technique.  Such  results  were 
analyzed  to  obtain  the  energy  differences  between  the  fully-on-glide-plane  core  structure, 
fully-on-cross-slip  plane  core  structure  and,  the  activated  structures  for  this  intersection. 

The  approach  used  to  obtain  activation  energies  for  transfonnation  from  one 
core  structure  to  another  is  as  follows:  As  the  stress  is  varied,  the  relative  magnitudes  of 
Escaig  stresses  on  the  glide  plane  and  cross-slip  plane  change,  thus  affecting  a  variation 
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in  the  length  of  the  partially  cross-slipped  region  of  the  screw-character  dislocation.  Note 
that  stabilizing  each  configuration  requires  manipulation  of  stresses  before,  after  and  at 
the  inflection  point  in  the  energy-distance  coordinate.  By  releasing  the  stress  from  the 
stabilized  state,  the  configuration  was  allowed  to  “fall”  towards  either  side  of  the  saddle 
point  configuration.  This  process  allowed  the  identification  of  the  saddle  point 
configuration  (within  the  energy  error  of  the  technique).  Once  identified,  the  absolute 
energy  of  this  configuration  without  the  applied  stress  is  calculated  by  freezing  a  small 
volume  of  atoms  near  the  core  transformation  region  and  then  minimizing  the  energy  of 
all  other  atoms  in  the  cell.  This  gives  an  upper  bound  for  the  energy  of  the  saddle  point. 
The  activation  energy  is  obtained  by  subtracting  the  energy  of  the  fully-glide-plane- 
spread  core  configuration  from  that  of  the  saddle  point  configuration.  The  primary  error 
in  this  calculation  arises  from  the  fact  that  an  arbitrary  number  of  atoms  are  kept  fixed 
while  obtaining  the  energy  of  the  saddle  point  configuration  without  applied  stress. 

2.1  Interatomic  Potential 

The  embedded  atom  potentials  used  for  the  simulations  are  the  potentials  developed 
for  FCC  Ni  by  Angelo,  Moody  and  Baskes  [8]  based  on  the  Voter  and  Chen  fonnat 
(Angelo),  as  well  as  the  two  Ni  potentials  developed  based  on  the  Voter  and  Chen  fonnat 
which  Rao  et.al.  used  in  their  bulk  cross-slip  simulations,  vni  and  vnih  [9].  Table  1  gives 
the  lattice  parameter,  cohesive  energy,  elastic  constants  and  stacking  fault  energy  for  each 
of  the  potentials.  The  three  Ni  potentials  used  in  the  simulations  give  almost  identical 
elastic  constants,  cohesive  energy  and  lattice  parameter,  whereas  the  stacking  fault 
energy  given  by  the  potentials  varies  from  60  -  120  mJ/rrr.  The  Shockley  partial  spacing 
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width,  d,  for  the  screw  dislocation  varies  from  d/b  =  5  for  potential  vnih  to  d/b  =8  for 
potential  vni.  The  Angelo  potential  gives  a  d/b  ratio  of  approximately  6.  This  allows  a 
detennination  of  cross-slip  activation  energy  as  a  function  of  d/b  ratio. 

2.2  Depiction  of  core  structures 

In  order  to  illustrate  the  relaxed  screw  dislocation  geometries  we  take  advantage  of 
the  increase  in  atomic  energy  produced  by  the  strain  field  of  the  partial  dislocations.  By 
plotting  the  atoms  having  assigned  energies  within  LAAMPs  of  greater  than  -4.42  or  - 
4.40  eV  (the  energy  of  atoms  in  the  stacking  fault  region)  the  partial  dislocations  can  be 
imaged  easily,  even  for  these  large  simulation  cells.  In  order  to  illustrate  the  cross- 
slipped-segment  products  of  the  screw  dislocation  the  positions  are  shown  in  a  [111]  as 
well  as  a  [112]  projection.  For  the  [112]  projection,  segments  spread  on  the  initial  (111) 
plane  appear  as  a  single  line  and  cross-slipped  segments  (i.e.  on  a  (1 1  1 )  plane)  appear  as 
a  pair  of  partials  separated  by  a  stacking  fault. 

3.0  Review  of  atomistic  core  structures  for  the  120°  intersection 

Figures  1-  3  depict  the  [111]  and  [112]  projections  (i.e.  the  „x-y’  and  „x-z’  planes)  of 
the  various  core  structures  obtained  for  the  screw  dislocation  -  120°  dislocation 
intersection,  discussed  in  reference  1,  obtained  using  the  Angelo  EAM  potential.  The 
starting  configurations  and  boundary  conditions  are  discussed  in  detail  in  the  literature 
[1].  Figure  1  illustrates  a  screw  dislocation  fully  on  the  glide  (111)  plane  and  constricted 
(i.e.  Stroh  constrictions  [10])  at  both  the  intersections  around  one  of  the  Shockley  partials 
of  the  120°  VZ[  1  01]  dislocation  residing  on  the  (111)  plane.  Figure  2  depicts  two  of  the 
core  structures  obtained  for  the  screw  dislocation  -120°  dislocation  intersection,  Lomer- 
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Cottrell  lock  configurations  (LC1  and  LC2)  in  the  [111]  and  [112]  projections.  After 
relaxation,  figures  2a  and  2b,  the  screw  dislocation  fully  resides  on  the  cross-slip  (111) 
plane.  In  figure  2b,  as  in  figure  1,  the  screw  dislocation  is  constricted  at  both  the 
intersections  around  one  of  the  Shockley  partials  of  the  120°  Vi[  101]  dislocation  residing 
on  the  (111)  plane.  In  figure  2a,  the  a/6[  1  12]  Shockley  partial  of  the  screw  dislocation 
combines  with  the  a/6[12  1  ]  Shockley  partial  of  the  120°  dislocation  to  fonn  an  extended 
node  resulting  in  a  short  segment  having  Burgers  vector  of  the  type  l/6[  0  1  1]  (e.g.  a 
Lomer-Cottrell  barrier).  At  the  other  end  of  the  Lomer-Cottrell  barrier,  the  screw 
dislocation  is  constricted  to  form  one  half  of  a  Stroh  constriction.  Figure  3  shows  the 
[111]  and  [112]  projections  for  two  of  the  core  structures  obtained  for  the  screw 
dislocation  -  120°  dislocation  intersection,  identified  as  Partially  Cross-Slipped  Locks 
(PCS1  and  PCS2).  In  both  cases,  the  screw  dislocation  is  partially  on  the  glide  (111) 
plane  and  partially  on  the  cross-slip  (111)  plane.  In  figure  3a,  as  in  figure  1,  the  screw 
dislocation  is  constricted  at  both  the  intersections  (negative  constrictions  [9,11])  around 
one  of  the  Shockley  partials  of  the  120°  Vi[  101]  dislocation  residing  on  the  (111)  plane. 
Since  the  core  is  in  the  partially  cross-slipped  state,  a  companion  positive  constriction  [9, 
11]  forms  between  the  intersections.  In  figure  3b,  as  in  figure  2b,  a/6[  1  12]  Shockley 
partial  of  the  screw  dislocation  combines  with  the  a/6[12  1  ]  Shockley  partial  of  the  120° 
dislocation  to  fonn  an  extended  node  resulting  in  a  short  segment  having  Burgers  vector 
of  the  type  l/6[  0  1  1]  (e.g.  a  Lomer-Cottrell  barrier).  At  the  other  end  of  the  Lomer- 
Cottrell  barrier,  the  screw  dislocation  is  constricted  to  fonn  one  half  of  a  negative 
constriction  and  is  cross-slipped  onto  the  glide  (111)  plane.  Since  the  core  is  in  the 
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partially  cross-slipped  state,  a  companion  full  positive  constriction  [9,  11]  forms  between 
the  intersections. 

3.1  Schematic  description  of  the  energy  well 

Figure  4  gives  a  schematic  plot  of  the  energy  of  the  screw  dislocation  configuration  at 
the  intersection  as  a  function  of  the  length  of  screw  dislocation  cross-slipped  onto  the 
cross-slip  plane  in  between  the  120°  intersecting  dislocations.  Previous  molecular  statics 
simulation  results  [1]  show  that  there  are  local  minima  in  the  energy  curve  when:  the 
screw  dislocation  length  on  cross-slip  plane  goes  to  zero  (fully-glide  state),  the  entire 
screw  dislocation  length  is  on  the  cross-slip  plane  (fully-cross-slip  state),  and  when 
approximately  half  the  screw  dislocation  length  between  the  intersections  is  on  the  cross¬ 
slip  plane  (partial-cross-slip  state).  Two  local  maxima  are  expected  for  the  energy  versus 
state  curve,  in  between  the  fully-glide  and  partial-cross-slip  state  and  the  fully-cross-slip 
and  partial-cross-slip  state.  Starting  from  the  partial-cross-slip  state,  application  of 
increasing  Escaig  stresses  on  the  configuration  are  used  to  drive  the  screw  dislocation  up 
the  energy  well  until  a  point  of  inflexion  (e.g.  the  first  derivative  of  the  energy  curve  and 
the  stress  is  a  maximum).  For  stresses  slightly  above  the  point  of  inflexion  in  the  energy 
curve,  the  screw  dislocation  will  collapse  to  either  the  fully-glide  or  fully-cross-slip  state. 
The  work  done  by  the  Escaig  stresses  up  to  the  inflexion  point  gives  the  energy  required 
for  the  screw  dislocation  to  transfer  from  the  partially-cross-slipped  state  to  the  fully- 
glide-plane  or  fully-cross-slip-plane  states.  The  maximum  Escaig  stresses  under  which 
the  partial-cross-slipped  state  is  stable  are  given  by  the  maximum  in  the  first  derivative  of 
this  energy  profile.  The  following  procedure  was  used  to  determine  the  local  maxima, 


that  is  the  activated  state  for  transferring  the  screw  dislocation  from  fully-glide-plane  or 
fully-cross-slip-plane  states  to  the  partially-cross-slipped  state  in  the  energy  curve. 
Partially-cross-slipped  states,  under  Escaig  stresses  slightly  above  the  maximum  in  first 
derivative  of  the  energy  curve,  are  relaxed  and  stopped  at  different  points  along  the 
energy  curve  and  the  Escaig  stresses  relaxed  back  to  zero.  Depending  on  the  initial  value 
of  the  Escaig  stress  and  the  position  on  the  activation  surface,  under  zero  stress  the  screw 
dislocation  configuration  will  collapse  to  one  of  the  three  possible  minimum  energy 
states:  partially-cross-slipped,  fully-glide-plane  or  fully-cross-slipped.  The  point  at 
which  the  screw  dislocation  transitions  from  relaxing  to  the  partially-cross-slipped  state 
to  relaxing  to  fully-glide-plane  or  fully-cross-slipped-plane  states  is  identified  as  the 
maxima  in  the  energy  curve  or  the  activated  state  for  cross-slip.  The  energy  of  this 
configuration  relative  to  the  fully-glide-plane  or  fully-cross-slip-plane  states  is  evaluated 
by  using  2  relaxation  techniques  from  the  activated  configuration,  which  are  described  in 
section  3.4  below. 

3.2  Behavior  of  cores  PCS1  and  PCS2  under  stresses 

Figure  5  gives  the  relaxed  core  structures  in  the  [112]  projection  for  the  screw  -  120° 
intersection  (core  structure  PCS1)  for  uniform  compressive  applied  stresses  of  250,  500, 
750,  1000,  1100  and  1150  MPa  along  the  [111]  direction,  obtained  using  the  Angelo 
EAM  potential.  For  the  1150  MPa  case,  an  additional  increment  of  50  MPa  compressive 
stress  was  applied  on  the  relaxed  1100  MPa  compressive  stress  configuration.  With 
increasing  compressive  applied  stress,  the  amount  of  cross-slipped  screw  segment  in 
between  the  two  intersecting  120°  dislocations  monotonically  increases  until  a  stress  of 
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1100  MPa  is  reached.  When  the  compressive  applied  stress  is  slightly  increased  from 
1 100  to  1150  MPa,  the  screw  dislocation  in  between  the  intersections  suddenly  collapses 
to  core  structure  LC2,  fully  residing  on  the  cross-slip  (111)  plane.  For  compressive 
applied  stresses,  the  effective  stacking-fault  energy  on  the  cross-slip  (111)  plane 
decreases  by  an  amount  equal  to  Tebe,  where  re  =  Sac,  S  is  the  Schmidt  factor  acting  on 
the  edge  component  of  the  Shockley  partials  on  the  cross-slip  (111)  plane,  ac  is  the 
applied  compressive  stress  along  the  [111]  direction  and  be  is  the  magnitude  of  the  edge 
component,  =  0.07 1  mn.  Similarly,  figure  6  gives  the  relaxed  core  structures  in  the  [112] 
projection  for  the  screw  -  120°  intersection  (core  structure  PCS1)  for  uniform  tensile 
applied  stresses  of  250,  500,  750,  1000  and  1050  MPa  along  the  [111]  direction,  again 
obtained  using  the  Angelo  EAM  potential.  As  before,  for  the  1050  MPa  case,  a  slight 
additional  increment  of  50  MPa  tensile  stress  was  applied  on  the  relaxed  1000  MPa 
tensile  stress  configuration.  In  this  case,  the  amount  of  cross-slipped  segment  in  between 
the  two  intersecting  120°  dislocations  montonically  decreases  until  a  stress  of  1000  MPa 
is  reached.  When  the  tensile  applied  stress  is  slightly  increased  from  1000  to  1050  MPa, 
the  screw  dislocation  in  between  the  intersections  suddenly  relaxes  to  core  structure  Gil, 
fully  residing  on  the  glide  (111)  plane.  For  tensile  applied  stresses,  the  effective 
stacking-fault  energy  on  the  cross-slip  (111)  plane  increases  by  an  amount  equal  to  xebe, 
where  tc  =  Sat,  at  is  the  applied  tensile  stress  along  the  [111]  direction.  Similar  behavior 
was  observed  for  core  structure  PCS2.  In  this  case,  when  the  applied  stress  was 
compressive,  core  structure  PCS2  relaxes  to  core  structure  LC1,  when  the  compressive 
stress  was  increased  from  650  to  700  MPa.  Similarly,  for  tensile  applied  stresses,  core 
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structure  PCS2  realxes  to  core  structure  Gil,  when  the  applied  tensile  stress  was 
increased  from  1 150  to  1200  MPa. 

3.3  Analysis  of  cores  PCS1  and  PCS2  under  applied  stresses 

Figures  7a  and  7b  illustrate  the  change  in  length  of  the  cross-slipped  segment  in 
between  the  two  intersecting  120°  dislocations  (x-axis)  versus  the  change  in  effective 
stacking-fault  energy  [y-axis]  due  to  the  applied  compressive  or  tensile  stresses  along  the 
[ill]  direction,  for  initial  core  structure  PCS1,  obtained  using  the  Angelo  EAM  potential. 
Figures  7a  and  7b  clearly  show  that  for  a  certain  critical  increase  or  decrease  in  stacking- 
fault  energy,  8y,  the  core  structure  for  the  intersection  collapses  to  either  LC2,  fully 
residing  on  the  cross-slip  plane  (compressive  stress)  or  Gil,  fully  residing  on  the  glide 
plane  (tensile  stress).  The  work  done  by  the  applied  stresses  up  to  the  critical  point  is 
given  by  the  area  under  the  curve  up  to  the  critical  point,  multiplied  by  the  Shockley 
partial  spacing  (estimated  to  be  6b,  where  b  is  the  magnitude  of  the  Burgers  vector  of  the 
screw  dislocation  [9])  and  is  equal  to  the  energy  difference  between  the  activated  states 
and  the  partially  cross-slipped  state,  PCS1,  8Eact-pcsi-  This  is  detennined  to  be  0.68  eV 
for  the  activated  state  on  the  compressive  side  and  1 .44  eV  for  the  activated  state  on  the 
tensile  side.  The  Escaig  stresses  required  for  transferring  the  dislocation  from  the 
partially-cross-slipped  state  (PCS1),  to  the  fully-glide-plane  state  (on  the  tensile  side, 
Gil),  is  314  MPa  and,  to  the  fully-cross-slipped  state  (LC2)  is  346  MPa. 

Similarly,  figures  8a  and  8b  show  the  change  in  length  of  the  cross-slipped  segment 
in  between  the  two  intersecting  120°  dislocations  (x-axis)  versus  the  change  in  effective 
stacking-fault  energy  [y-axis]  due  to  the  applied  compressive  or  tensile  stresses  along  the 
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[Ill]  direction,  for  initial  core  structure  PCS2,  obtained  using  the  Angelo  EAM  potential. 
The  Escaig  stresses  required  for  transferring  the  dislocation  from  the  partially  cross- 
slipped  state  (PCS2)  to  the  fully-glide-plane  state  (on  the  tensile  side,  GL1),  is  361  MPa 
and  to  the  fully-cross-slipped  state  (LC1)  is  204  MPa. 

3.4  Determination  of  cross-slip  activation  energy 

To  determine  the  activation  energy  for  cross-slip  from  the  fully-glide-plane  state  (Gil) 
to  the  partially-cross-slipped  state  (PCS1),  the  simulation  having  1050  MPa  tensile  stress 
applied  on  the  partially-cross-slipped  state  PCS1,  was  stopped  after  varying  amounts  of 
relaxation  and  the  applied  stress  relaxed  back  to  zero.  Conjugate  gradient  minimization 
after  the  stress  is  relaxed  back  to  zero  indicated  whether  the  core  relaxes  to  the  partially- 
cross-slipped  state  (PCS1)  or  the  fully-glide-plane  state  (Gil).  The  configuration  at 
which  the  relaxation  transitions  from  relaxing  to  PCS1  to  relaxing  to  Gil,  is  taken  to  be 
the  activated  state  for  cross-slip  from  Gil  to  PCS1.  Figure  9  gives  the  structure  of  the 
activated  state  in  the  [112]  projection  for  the  screw  -  120°  intersection  (core  structure 
PCS1),  obtained  using  the  Angelo  potential,  which  gives  a  stacking-fault  energy  of  90 
mJ/m  .  Figure  9a  shows  that  in  the  activated  configuration,  a  small  portion  of  the  screw 
dislocation  near  the  intersection,  of  length  ~1.5  d,  resides  on  the  cross-slip  plane.  The 
attractive  interaction  between  the  negative  and  positive  constrictions  goes  to  zero  at  a 
separation  distance  of  1.5  d,  beyond  which  it  becomes  repulsive,  in  the  presence  of  an 
intersecting  dislocation.  This  behavior  is  significantly  different  from  typical  Escaig 
process  for  cross  slip  where  the  attractive  interaction  between  the  constrictions  only  goes 
to  zero  at  an  infinite  separation  distance  [3-5],  Similar  analysis  was  perfonned  to 
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determine  the  activation  energy  for  cross  slip  from  the  fully-cross-slip-plane  state  (LC2) 
to  the  partially-cross-slipped  state  (PCS1).  In  this  case  the  simulation  having  1150  MPa 
compressive  stress  applied  on  the  partially  cross-slipped  state  PCS1  was  stopped  after 
varying  amounts  of  relaxation  and  the  applied  stress  relaxed  back  to  zero.  Conjugate 
gradient  minimization  after  the  stress  is  relaxed  back  to  zero  indicated  whether  the  core 
relaxes  to  the  partially-cross-slipped  state  (PCS1)  or  the  fully-cross-slip-plane  state 
(LC2).  The  configuration  at  which  the  relaxation  transitions  from  relaxing  to  PCS1,  to 
collapsing  to  LC2,  is  taken  to  be  the  activated  state  for  cross-slip  from  LC2  to  PCS1. 
Figure  9b  shows  that  in  the  activated  configuration,  a  small  portion  of  the  screw 
dislocation  near  the  intersection,  of  length  ~d,  resides  on  the  glide  plane.  The  attractive 
interaction  between  the  negative  and  positive  constrictions  goes  to  zero  at  a  separation 
distance  of  d  beyond  which  it  becomes  repulsive,  in  the  presence  of  an  intersecting 
dislocation.  Such  a  behavior  is  significantly  different  from  the  standard  Escaig  process. 

The  activation  energy  for  cross-slip  from  state  Gil  to  state  PCS1,  8EGu-pcsi  was 
detennined  using  two  different  relaxation  techniques  from  the  identified  activated 
configuration.  In  the  first  technique,  the  energy  was  detennined  for  relaxation  from  the 
activated  state,  to  state  Gil,  under  an  applied  tensile  stress  of  1050  MPa  along  the  [111] 
direction.  This  energy  was  corrected  for  two  factors  to  obtain  the  activation  energy  for 
cross-slip:  a)  work  done  by  the  applied  Escaig  stresses  in  collapsing  the  activated 
configuration  (of  the  order  of  0.2  eV)  and,  b  the  increase  in  energy  of  the  activated 
configuration  due  the  applied  Escaig  stresses  (of  the  order  of  0.4  eV).  For  the  second 
technique,  the  activated  state  was  relaxed  to  state  Gil  having  zero  applied  stress,  with  and 
without  atoms  within  a  cube  of  size  1.2  nm  centered  at  the  positive  constriction  of  the 
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activated  state  fixed.  The  difference  in  relaxation  energies  between  the  two  types  of 
simulations  was  taken  to  be  the  activation  energy  for  cross  slip  from  Gil  to  PCS1.  The 
results  for  the  cross-slip  activation  energy  obtained  using  these  two  techniques  of 
relaxation  were  within  0.2  eV  of  each  other  and,  are  taken  to  be  representative  of  the 
error  in  the  determination  of  the  activation  energy  for  cross  slip  using  this  procedure. 
Similar  simulations  and  analysis  were  perfonned  to  detennine  the  activation  energy  for 
cross-slip  from  Gil  to  PCS1  using  two  other  Ni  EAM  potentials,  vni  and  vnih.  Also, 
similar  simulations  and  analysis  were  performed  to  detennine  the  activation  energy  for 
cross  slip  from  the  fully-cross-slip-plane  state  (LC2)  to  the  partially-cross-slipped  state 
(PCS1)  using  all  three  EAM  potentials,  8ELC2-pcsi-  In  this  case,  the  two  techniques  of 
relaxation  described  before  for  detennining  the  cross-slip  activation  energy  gave  values 
that  differed  by  less  than  0.1  eV.  Figure  10  is  a  plot  of  the  cross-slip  activation  energies 
for  cross  slip  from  state  Gil  to  state  PCS1,  as  well  as  from  state  LC2  to  state  PCS1, 
obtained  using  the  three  Ni  EAM  potentials,  versus  the  parameter  (d/b)ln(d/b)  given  by 
the  potentials.  Also  plotted  in  figure  10  are  the  corresponding  values  for  the  cross-slip 
activation  energy  determined  from  simulations  of  the  Escaig  process  [ref].  Figure  10 
clearly  shows  that  the  cross-slip  activation  energy  is  linearly  proportional  to  (d/b)ln(d/b) 
both  at  the  intersection  and  for  the  Escaig-type  process.  The  cross-slip  activation-energy 
(from  Gil  to  PCS1)  values  near  the  intersection  are  a  factor  of  3-4  lower  than  that  for  the 
Escaig  process.  Also,  the  cross-slip  activation  energy  from  the  fully-cross-slip-plane 
state  (LC2)  to  the  partially-cross-slipped  state  (PCS1),  is  approximately  a  factor  of  6 
lower  than  that  for  the  Escaig  process.  The  8ELc2-pcsi  values  can  be  identified,  by 
symmetry,  as  the  cross-slip  activation  energy  from  the  fully-glide  plane-state  to  the 
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partially-cross-slipped  state  at  a  120°  screw  dislocation  intersection,  where  the 
intersecting  forest  dislocation  has  a  1/2<011>  Burgers  vector,  <101  >  line  direction  and 
resides  on  the  (111)  plane.  This  suggests  that  the  cross-slip  activation  energy  near  an 
intersection  is  strongly  dependent  on  the  Burgers  vector  of  the  intersecting  dislocation 
and  could  be  a  factor  of  6  lower  than  that  usually  considered.  To  obtain  the  cross-slip 
activation  energies  dislocations  intersecting  forest  dislocations  in  other  FCC  metals  from 
Fig.  10,  one  has  to  scale  by  the  elastic  constants  of  the  metal  relative  to  that  of  Ni  at  the 
appropriate  d/b  ratio  for  that  metal.  Note  that  Fig.  10  shows  that  intersection  cross-slip 
activation  energies  for  typical  FCC  metals  such  as  Ni,  Cu  and  Au  ranges  from  0.3  -  0.7 
eV.  The  activation  volume  for  the  cross-slip  process  with  respect  to  applied  Escaig 
stresses  is  given  by  the  area  of  the  cross-slipped  portion  of  the  screw  character  dislocation 
in  the  activated  state  (Fig.  9)  multiplied  by  the  magnitude  of  the  edge  component  of  the 
Burger’s  vector  of  the  Shockley  partials.  From  the  simulations,  values  for  this  parameter 
range  from  10  -  20b  for  cross-slip  at  screw  dislocation  intersections. 

4.0  Discussion 

The  results  of  this  work  clearly  show  that  cross-slip  nucleation  at  screw 
dislocations  intersecting  forest  dislocations  is  relatively  easy,  having  an  activation  energy 
that  is  a  factor  of  3-6  lower  than  that  for  the  Escaig  process.  This  is  an  important  result 
since  the  probability  of  dislocation  intersections  is  relatively  large  as  compared  to  other 
mechanisms  of  obstacle  induced  cross-slip,  like  cross-slip  at  screw  dipoles,  jogs,  and 
surfaces  [13-15].  Further,  cross-slip  at  screw  dipoles,  jogs  and  surfaces  is  athermal  only 
under  specific  conditions  [13-15],  The  current  result  rationalizes  the  profuse  nature  of 
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cross  slip  in  FCC  crystals,  unlike  the  current  models  of  thermally-activated  cross  slip  that 
require  high  stress  to  provide  a  self-consistent  explanation.  This  finding  should  allow 
higher-level  mesoscale  models  of  dislocation  behavior  to  better  represent  the  cross-slip 
process  without  resort  to  ad  hoc  assumptions  about  obstacles.  It  should  be  noted  that 
since  a  molecular  statics  technique  has  been  used  in  this  work,  the  cross-slip  activation- 
energy  values  reported  in  this  manuscript  are  expected  to  be  upper  bounds.  A  more 
sophisticated  molecular  dynamics  reaction  pathway  technique  [16]  might  provide  more 
accurate  values  for  the  activation  energies  of  cross  slip  near  forest  dislocation 
intersections. 

The  new  cross-slip  nucleation  model  has  a  variety  of  implications  for  crystal 
plasticity  in  FCC  materials.  For  example,  within  the  present  mechanism,  the  frequency 
of  cross-slip  should  scale  with  the  forest  dislocation  density.  The  growth  of  such  nuclei 
should  depend  upon  the  relative  magnitude  of  local  stresses  on  the  glide  plane  and  the 
cross-slip  plane  at  the  partially  cross-slipped  screw-dislocation  intersection  region.  Such 
behavior  of  the  partially  cross-slipped  core  under  different  modes  of  applied  stress  could 
be  studied  using  atomistic  simulations  as  well  as  dislocation  dynamics  simulations. 
Similar  activation  analysis  using  atomistic  simulations  must  be  performed  for  other 
intersections  (varying  line  direction  of  the  intersecting  dislocation)  on  the  (111)  plane 
that  form  Glide,  Lomer-Cottrell  or  Hirth  locks  to  determine  the  dependence  of  the 
activation  energy  for  cross-slip  on  the  line  direction  as  well  as  the  Burger’s  vector  of  the 
intersecting  forest  dislocation.  Also,  it  would  be  instructive  to  perform  activation 
analysis  using  atomistic  simulations  of  these  screw-dislocation  intersection  core 
structures  for  other  FCC  materials  having  different  stacking  fault  energies  (i.e.  A1  and 
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Cu)  to  verify  the  results  in  Ni.  The  intersection  mechanism  of  cross-slip  nucleation 
should  also  be  implemented  in  3d  dislocation  dynamics  simulations  as  an  alternative  to 
Escaig’s  model  for  FCC  materials.  Finally,  it  would  be  instructive  to  revisit  the 
Bonneville  and  Escaig  experimental  results  on  cross-slip  in  FCC  Cu  in  light  of  the  new 
intersection  mechanism  for  cross-slip  nucleation  [2,  17]. 

According  to  Washburn  [18],  double-intersection  cross  slip,  where  the  segment 
that  has  been  pulled  into  the  cross-slip  plane  soon  encounters  another  attractive 
intersection  that  brings  it  back  onto  another  primary  glide  plane,  provides  a  reasonable 
mechanism  for  dislocation  multiplication  and  the  growth  of  slip  bands  at  low 
temperatures.  Naturally,  such  a  mechanistic  process  must  be  demonstrated  via  modem 
simulation  methods.  Also,  classical  theories  of  strain-hardening  assume  that  dislocation 
storage  in  Stage  II  of  single-crystal  FCC  materials  is  a  result  of  junction  formation  [19] 
or,  2-dimensional  concave  loop  formation  [20]  as  the  gliding  dislocation  traverses 
through  an  array  of  forest  dislocation  obstacles  on  its  glide  plane.  However,  one  of  the 
major  problems  in  classical  strain-hardening  models  is  to  explain  how  the  generation  of  a 
three-dimensional  network  of  stored  dislocations  occurs  as  a  consequence  of  two- 
dimensional  glide  [20].  We  note  that  the  intersection  cross-slip  nucleation  mechanism 
for  dislocation  storage  may  provide  a  convenient  mechanism  of  generating  a  three- 
dimensional  network  of  stored  dislocations  from  two-dimensional  glide. 

Finally,  the  anomalous  flow  behavior  of  intermetallics  (i.e.  M3AI)  has  thus  far  been 
modeled  using  core  transfonnations  that  require  thennal  activation.  Here  again,  the  key 
issue  has  been  that  constriction  energies  are  prohibitively  high  [3,4,21],  The  possibility 
of  the  core  transformation  in  the  presence  of  interaction  with  other  dislocations  provides 
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a  novel  mechanism  to  consider  and  explore  in  rationalizing  the  flow  anomaly  of  these 
types  of  materials. 

5.0  Summary 

Atomistic  simulations  were  conducted  to  identify  the  saddle-point  configuration 
and  thus  the  activation  energy  for  cross  slip  of  a  screw  dislocation  in  the  presence  of  a 
pair  of  forest  dislocations.  The  approach  involved  the  use  of  external  stress  along  [ill] 
to  vary  the  segment  length  for  a  partially  cross  slipped  VifllO]  screw  dislocation 
intersecting  a  pair  of  120°  Vi[  101]  Burgers  vector,  <01  1  >  line  direction  dislocations 
residing  on  the  (111)  plane.  The  simulations  in  this  work  show  that  the  nucleation  of 
cross-slip  at  dislocation-dislocation  interactions  is  significantly  more  probable  than  the 
Escaig  process,  having  an  activation  barrier  that  is  a  factor  of  3  -  6  lower.  Further,  this 
work  shows  that  atomistic  simulations  can  be  used  to  identify  types  of  interactions  that 
result  in  cross-slip  and  these  results  in  turn  may  be  used  in  higher-level  mesoscale 
simulations.  Summarizing  the  results:  equilibrium  partially  cross-slipped  core  structures 
under  applied  compressive  or  tensile  stresses  along  the  [111]  direction  were  determined 
using  atomistic  simulations.  Beyond  a  critical  stress,  the  equilibrium  configurations  relax 
to  a  core  structure  that  is  either  fully  on  the  glide  plane  for  tensile  stresses  or  fully  on  the 
cross-slip  plane  for  compressive  stresses. 

1)  The  activation  barrier  for  cross-slip  from  the  fully-glide-plane  state  (Gil)  to  the 
partially-cross-slipped  state  (PCS1),  8Egii-pcsi  at  these  intersections,  is  determined  to  be  a 
factor  of  3-4  lower  than  that  for  the  Escaig  process.  While  the  new  cross  slip  mechanism 
has  a  lower  activation  energy  it  has  the  same  linear  relationship  to  (d/b)ln(d/b)  as  the 


18 


Escaig  process.  The  activation  barrier  for  cross-slip  from  the  fully-cross-slip-plane  state 
(LC2)  to  the  partially-cross-slipped  state  (PCS1),  8ELC2-pcsi  at  these  intersections,  is 
detennined  to  be  approximately  a  factor  of  6  lower  than  for  the  Escaig  process.  These 
8Elc2-pcsi  values,  by  symmetry,  are  identified  as  the  activation  barrier  for  cross-slip  from 
the  fully- glide-plane  state  to  the  partially-cross-slipped  state  at  a  120°  screw  dislocation 
intersection,  where  the  intersecting  dislocation  has  a  'A<0 1 1  >  Burgers  vector,  <  10  1  >  line 
direction  and  resides  on  the  (111)  plane. 

2)  Cross-slip  activation-energy  values  near  forest  dislocation  intersections  in 
typical  FCC  metals  like  Ni,  Cu  and  Au  range  from  0.3  -  0.7  eV. 
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Figure  captions 


Figure  1 :  [111]  (X-Y  plane)  and  [112]  (X-Z  plane)  projections  of  the  core  structure  (Gil) 
for  a  screw- 120°  intersection  in  FCC  Ni  from  atomistic  simulations.  Atoms  having  an 
energy  greater  than  the  energy  at  a  stacking  fault  are  plotted.  The  axes  dimensions  are  in 
units  of  Angstroms. 

Figure  2:  [111]  (X-Y  plane)  and  [112]  (X-Z  plane)  projections  of  the  core  structures 
(LC1  and  LC2)  for  a  screw-1200  intersection  in  FCC  Ni  from  atomistic  simulations. 
Atoms  having  energies  greater  than  the  energy  at  a  stacking  fault  are  plotted.  The  axes 
dimensions  are  in  units  of  Angstroms. 

Figure  3:  [111]  (X-Y  plane)  and  [112]  (X-Z  plane)  projections  of  the  core  structures 
(PCS1  and  PCS2)  for  a  screw- 120°  intersection  in  FCC  Ni  from  atomistic  simulations. 
Atoms  having  energies  greater  than  the  energy  at  a  stacking  fault  are  plotted.  The  axes 
dimensions  are  in  units  of  Angstroms. 

Figure  4:  a)  A  schematic  plot  of  the  energy  profile  of  the  screw-dislocation  configuration 
as  a  function  of  the  length  of  the  screw  dislocation  on  the  cross-slip  plane  in  between  the 
forest  dislocation  intersections,  b)  A  schematic  sketch  of  the  [112]  projection  of  the 
screw-dislocation  configuration  in  between  forest  dislocation  intersections  under  applied 
Escaig  stresses. 
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Figure  5:  The  [112  ]  projection  of  the  equilibrium  configurations  of  the  partially  cross- 
slipped  core  structure,  PCS1,  under  uniform  applied  compressive  stresses  of  0,  250,  500, 
750,  1000,  1 100  and  1150  MPa  applied  along  the  [111]  direction.  Atoms  having  energies 
greater  than  the  energy  at  a  stacking  fault  are  plotted.  The  axes  dimensions  are  in  units  of 
Angstroms. 

Figure  6:  The  [112]  projection  of  the  equilibrium  configurations  of  the  partially-cross- 
slipped  core  structure,  PCS1,  under  uniform  applied  tensile  stresses  of  0,  250,  500,  750, 
lOOOand  1050  MPa  applied  along  the  [111]  direction.  Atoms  having  energies  greater 
than  the  energy  at  a  stacking  fault  are  plotted.  The  axes  dimensions  are  in  units  of 
Angstroms. 

Figure  7:  A  plot  of  the  increase  or  decrease  in  length  of  the  cross-slipped  segment  in 
between  the  two  intersecting  120°  dislocations  (x-axis)  versus  the  decrease  or  increase  in 
effective  stacking-fault  energy  [y-axis]  due  to  the  applied  compressive  or  tensile  stresses 
along  the  [111]  direction,  for  initial  core  structure  PCS1.  The  activation  barrier  for  the 
transfer  of  the  dislocation  from  the  partially-cross-slipped  state  (PCS1)  to  the  either  the 
fully-cross-slip  plane  state  (LC2)  or  the  fully-glide  plane  state  (Gil)  is  also  indicated. 

Figure  8:  A  plot  of  the  increase  or  decrease  in  length  of  the  cross-slipped  segment  in 
between  the  two  intersecting  120°  dislocations  (x-axis)  versus  the  decrease  or  increase  in 
effective  stacking-fault  energy  [y-axis]  due  to  the  applied  compressive  or  tensile  stresses 
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along  the  [111]  direction,  for  initial  core  structure  PCS2.  The  activation  barrier  for  the 
transfer  of  the  dislocation  from  the  partially-cross-slipped  state  (PCS2)  to  the  either  the 
fully-cross-slip  plane  state  (LC1)  or  the  fully-glide  plane  state  (Gil)  is  also  indicated. 

Figure  9:  The  [112  ]  projection  of  the  structure  of  the  activated  configuration  for  transfer 
of  screw  dislocation  from  a)  state  Gil  to  PCS1  (tensile  side)  and  b)  state  LC2  to  state 
PCS1  (compressive  side)  obtained  using  the  Angelo  EAM  potential  for  Ni.  The  x  axis  is 
along  the  [110]  direction  and  the  y-axis  is  along  the  [ill]  direction.  The  axes 
dimensions  are  in  units  of  Angstroms. 

Figure  10:  A  plot  of  the  activation  energy  for  cross-slip  from  state  Gil  to  state  PCS1 
(8Egii-pcsi),  as  well  as  from  state  LC2  to  state  PCS1  (8ELc2-pcsi),  in  Ni  as  a  function  of 
the  parameter  d/bln(d/b)  given  by  the  EAM  potentials,  Angelo,  vni  and  vnih.  Also  shown 
are  the  corresponding  cross-slip  activation  energy  values  for  the  Escaig  process. 
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Property 

Value  (Angelo) 

Value  (vni) 

Value  (vnih) 

a0 

0.352nm 

0.35180789nm 

0.35256848nm 

Ci, 

2.46  x  10"  N/m- 

2.44  x  1011  N/m" 

2.44  x  10'  ‘N/m2 

C12 

1.47  x  10“  N/m2 

1.49  x  10"  N/m2 

1.49  x  10H  N/m2 

C44 

1.25  x  10“  N/m2 

1.25  x  10“  N/m" 

1.19  x  10“  N/m2 

Ec 

-4.45  eV 

-4.43  eV 

-4.43  eV 

Y 

.089  J/nr 

.059  J/m" 

.119  J/m- 

Table  1:  Lattice  parameter,  a0;  Cohesive  energy,  Ec;  Elastic  constants  Cn,  C12  and  C44; 
and  stacking  fault  energy,  y,  given  by  the  3  Ni  EAM  potentials,  Angelo,  vni  and  vnih. 
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